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The most pecuUar radio characteristics of the TeV emitting high-mass X-ray binary 
LS I -|-61°303 are two periodicities: A large periodic outburst which exhibits the same period 
as the orbit (phase "!>) and a second periodicity of 1667 days (phase Q) which modulates the orbital 
phase and amplitude of the large outburst. Recent analysis of the radio spectral index present strong 
evidence for the presence of the critical transition from optically thick emission (related to a steady 
jet) to an optically thin outburst (related to a transient jet) as in other microquasars. In parallel, 
a switch from a low/hard X-ray state to a transitional state (e.g. steep power law state) would be 
expected. We show how the critical transition from optically thick emission to an optically thin 
outburst is modulated by O. Folding over a too large Q interval mixes up important information 
about the outbursts and can yield a false picture of the emission behaviour of the source along the 
orbit. We therefore analyse the implications of this long period for treatment of hard X-ray/high 
energy data obtained from LS I -|-61°303 , e.g. with Fermi-LAT or INTEGRAL, taking into account 
this long-term periodicity. 



I. INTRODUCTION 

LS I -|-61°303 is an X-ray binary formed by a com- 
pact object and a massive star with an optical spec- 
trum typical for a rapidly rotating BO V star [l| . The 
nature of the compact object, if neutron star or black 
hole, is still unknown due to the large uncertainty in 
the inclination of the object 0, It travels around 
its companion star on an eccentric orbit with a period 
of 26.5 days [4]. 

Radio spectral index analysis by Q have found two 
peaks along the orbit of LS I -|-61°303. Each peak 
shows the microquasar characteristic of a switch from 
a steady (optically thick) to a transient (optically 
thin) jet. Furthermore, each peak in the radio spectral 
index is accompanied by two distinguishable peaks in 
the radio flux. This confirms that there are really two 
different outbursts (an optically thick and an optically 
thin one) . Also high energy observations with EGRET 
[H and Fermi-LAT Q indicate two peaks along the 
orbit. The high energy peaks are supposed to be due 
to inverse Compton upscattering of the UV photons 
from the donor star by the relativistic electrons of the 
jet, which in turn strongly attenuates the radio peak 
around periastron (see discussion in [1] and 0). 

It has been shown from theory that for a micro- 
quasar with an eccentric orbit (LS I +61°303 : e=0.54- 
0.7 [1, Q), indeed, the different relationship between 
the accretion rate for density and velocity described 
by p^ . creates two peaks in the accretion rate curve, 
one at peri astron and a second one towards apastron 

The clear periodicity in LS I -|-61°303 is amongst its 
most peculiar characteristics. The orbital period ($) 
modulates the flux at radio, H-alpha, X-rays and also 
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FIG. 1: Radio light curve of LS I -h61°303 from the Green 
Bank interferometer vs. orbital phase 0]. The large out- 
burst towards apastron is clearly observable. The broad- 
ness of the curve stems from the 4.6 yr period which mod- 
ulates the orbital occurrence of the peak. 



in gamma-rays (high and very high energy) P. [gI [tI. [isl. 
|l6, 17, 18, 19, 20, 21]. But LS I -K61°303 holds another 
peculiarity: A second period, which modulates the 
radio flux over a long period of 6=4.6 yr 0]. This 
modulation in radio is shown in Fig. [1] (radio flux 
density vs. $), which shows how the peak flux shifts 
around apastron, and in Fig 21 where the radio flux 
density is given vs. Q. In addition, this long period 
has also been shown in H-alpha (see Fig. 2 in [22j). 
It is suggested that this modulation could be due to 
periodic she fl ej ections from the circumstellar disk of 
the Be star 1231. 
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FIG. 2: Top: Radio spectral index data for 6=0.788-0.927 
from GBI [y| . Bottom: _Fermi-LAT lightcurve obtained by 
0]. The two vertical lines indicate the peaks around pe- 
riastron and around apoastron. They both correspond to 
a negative spectral index in the radio data, which corre- 
sponds to optically thin emission in both cases. 
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FIG. 3: Top: Radio spectral index data for 0=0.6-0.85 
from GBI ^9j. Bottom: Radio spectral index data for 
0=0.1-0.3 from GBI ,2al. 



II. HIGH ENERGY OBSERVATIONS AND 
THE RADIO SPECTRAL INDEX 

The radio spectral index analysis by [Hj have shown 
the importance of the large phase 9 for the analysis of 
radio data from LS I -|-61°303 . It has thereby proven 
that the radio outbursts really consist of two consec- 
utive outbursts. These outbursts have completely dif- 
ferent characteristics. In the microquasar model, an 
optically thick outburst in radio is associated with a 
steady jet, centered on the compact object. The op- 
tically thin outburst comes from a transient jet, de- 
tached from the central engine. In the unified model 
of X-ray states with radio jets, the radio states in 
microquasars are clearly associated with two X-ray 
states. A steady jet corresponds to the low hard state, 
where the X-ray spectrum is characterized by T wl.5, 
whereas a transient jet corresponds to a transitional 
state, e.g. the steep power law state with F >2.4 
[24l I25I l2a | . Furthermore, high energy and very high 
energy emission are directly connected to the tran- 
sitional state (steep power law state) as the power 
law is without cut-off and extends into the gamma-ray 
regime. As a matter of fact, when LS I -|-61°303 was 
detected by MAGIC and VERITAS the spectrum was 
always fitted with a power law with an index F >2.4 
[13, [23, [13 ■ Moreover, as discussed in the next 
section, INTEGRAL observations seem to indicate a 



change of the photon index consistent with a change 
from the low hard to the transitional state in agree- 
ment with the radio spectral index. 

If the radio spectral index now tells us about the na- 
ture of the outburst, then the high and very high en- 
ergy observations should corroborate this nature and 
might give additional information about the emission 
processes. But in order to get the complete informa- 
tion from the radio spectral index for LS I -1-61° 303 , 
it was neccessary to take into account the long period 
Q. Therefore, to compare high energy data with the 
radio spectral index, it is neccessary as well to take 
only data of the same 9 phase. Q has done so e.g. 
for the first eight months of Fermi-LAT observations 
and the results (see Fig. [2]) show that in this com- 
parison the high energy peaks established by Fermi- 
LAT both correspond to the optically thin outbursts 
of LS I -|-61°303. One could then draw the conclu- 
sion that the population of relativistic particles, pro- 
ducing the optically thin outburst (and therefore the 
transient jet), are also responsible for the production 
of the high energy emission. Furthermore, in Fig. [3] 
are shown radio spectral index curves vs. the orbital 
phase for two different 9 intervals. In particular the 
orbital occurrence of the second peak is different in 
the different intervals. It is evident that folding data 
only on the orbital period, without respect to the long 
period, might result in mixing up different ejection 
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FIG. 4: Published high energy observations of 
LS I +61°303 by EGRET, MAGIC, fermi-LAT, VERI- 
TAS and INTEGRAL shown in the context of the 4.6 yr 
radio period (O) [29| . The underlying light curve gives 
the radio flux at 8.3 GHz against Q [5|. The minimum 
lies around O ~0.5 and the maximum around ~0.97. The 
three largest O intervals where the source was observed 
with INTEGRAL are named here II, 12 and 13. 



processes, because the orbital phase of the peaks and 
therefore of the switch from optically thick to thin 
emission does not stay the same over the 4.6 yr pe- 
riod. This point will be discussed in the next section. 



III. INTEGRAL OBSERVATIONS OF 
LS I +61°303 AND 6 

In Fig. 131 high energy and very high energy ob- 
servations of LS I +61°303 with different instruments 
are shown with respect to the Q intervals in which 
they have been carried out. Both, MAGIC and VER- 
ITAS probe the minimum and the maximum parts 
of [23, US] • Due to the long integration times 
for the data, the sampling within these intervals is 
not very strong, although EGRET observed one com- 
plete orbital cycle at each of the two given intervals 

HO, HH . The latest published coverage by Fermi is 
of 2.5 years fsO^. 

For INTEGRAL, the long-term monitoring of 
LS I +61°303 covers vast parts of the 4.6 yr period. 
We have denoted the three big INTEGRAL observa- 
tions with respect to 9 by II, 12 and 13, as noted in 
Fig. m Folding the data by orbital phase, of course, 
increases the sampling. By doing so, it was estab- 
lished that the emission between 10-100 keV is clearly 
modulated by the orbital phase [3, 0, HH . 

Following the radio spectral index though (see Fig. 
|3]) , folding over almost a complete Q cycle would imply 
that even though the resulting light curves could show 
the overall periodicity of the source at these energies, a 



deeper insight about the ejections can be obscured by 
mixing different ejection processes (namely optically 
thick and thin ejections). The spectral analyses can 
then get corrupted, because of the same mixing of 
ejection processes. 

II and 13 are covering most of the maximum of 9 
(see Fig. 2]). An interpretation of these data with 
respect to the flux and the spectrum should be less 
corrupted than a mix of data from the maximum and 
the minimum. In fact, [l8l |. who used most of 13 and 
parts of II, have found not only the modulation of the 
lightcurve with the orbital period, but also found that 
along the orbit the spectral index F changed from «1.5 
around periastron to «3.2 around apastron. This re- 
sult strengthens the two-peak accretion model, as for 
the 0-intervals covered by these observations an opti- 
cally thin ejection (a transient jet) is expected around 
apastron (see Q and discussions therein) As men- 
tioned above, in the unified X-ray states model with 
radio jets, a transient is associated with a transitional 
state, e.g. the steep power law, characterized by a 
spectral index F >2.4 [H HI, [2i|. 



IV. CONCLUSIONS AND DISCUSSION 

In the high mass X-ray binary LS I -|-61°303, two 
clear radio periodicities are present, one coincident 
with its orbital period (see Fig. [1]) and the other mod- 
ulating the strength of the large radio outburst over a 
period of 4.6 yr [J, [l^j. Both periods have also been 
observed in H-alpha emission. This system is amongst 
a few to have been detected not only in radio and X- 
rays, but also at high and very high energies. It is from 
these observations that, together with radio spectral 
index analysis, important insights into the nature of 
the system and its emission processes can be deduced. 
The radio spectral index tells us about the nature of 
the observed outburst. There are two different kinds 
of outburst: optically thick (spectral index a > 0) and 
optically thin (a < 0). The first is attributed in the 
microquasar model to a slow outfiow, while the latter 
is associated with an ultrarelativistic transient jet. 

Recent spectral index analysis of LS I -|-61°303 show 
that a subsequent realization of these two types 
of outburst takes places not only once, but twice 
along the orbit (see Fig. 2 top). This can be ex- 
plained by the two peak accretion/ejection micro- 
quasar model for systems with an eccentric orbit like 
LS I +61°303 (e=0. 54-0.7) (:5j and references within). 
In this model an expected radio peak around perias- 
tron is attenuated by inverse Compton losses due to 
the dense stellar UV field. In the first eight months 
observations with Fermi-LAT, a related peak at high 
energies is observed (see Fig. 2 at $=0.3-0.4). A 
second (but smaller) high energy peak is observed co- 
incident with the large radio outburst around apas- 
tron ($=0.65-0.75). Furthermore, as discussed in Q, 
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the two peak shape of the a vs. $ curve varies with 
0, as does the distance of the two peaks (see also 
Fig. 2 here). The two-peak microquasar model of [T^ 
predicts these variations in the accretion curve by in- 
corporating changing wind velocities for the Be star 
(see their Fig. 6). This behaviour should be seen at 
other wavelengths as well, due to the connected emis- 
sion processes. In fact, after the first eight months, 
-Fermz-LAT detected an increase in the overall flux 
level (O « 0.92) and a broadening of the peak shape 
[30| . These variations are consistent with the obser- 
vation that a varies with O. Variations are also seen 
with VERITAS and MAGIC. LS I -l-61°303 was de- 
tected around apastron until 2008 and then became 
quiescent between 2008-2010, where no detection was 
reported by VERITAS, while MAGIC reported only 
weak detection around apastron. Then in October 
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2010, VERITAS detected the source again, this time 
around periastron pTI. |27|. [28} . 

The strong connection between radio and HE/VHE 
emission underlines therefore the importance of the 
radio periods for the analysis of e.g. INTEGRAL data 
and can and should in principle be extended to other 
high energy instruments. 
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